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In spite of the fact that the two-body B decays into Xc 2 such as B ^ 
Xc2K^*^ are suppressed by the QCD factorization effect, the inclusive B —)■ 
Xc 2 X branching fraction amounts to one third of the non-suppressed B —)■ 
XciX decays because of the decay modes to the multi-body final states. 
Using a large statistics T(4S') data sample corresponding to 772 million B 
meson pairs accumulated by the Belle detector at the KEKB e’''e“ collider, 
precise measurements of inclusive B — )■ Xci and Xc2 branching fractions 
are carried out. The multi-body hnal states such as XcjKt^i XcjKt^t^ and 
so on are also investigated to look for new charmonium-like resonance. 
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1 Introduction 


Inclusive production of Xc2 mesons in B decays is relatively large Bia in spite of 
the fact that two-body B decays into Xc 2 are highly suppressed [31 HI [5] (due to 
the angular momentum conservation). Differential branching fraction {DB) in bins 
of Xcj {J = 1,2) [ 6 ] suggests that Xc 2 is found to be coming from three-body or 
higher multiplicity decays inEi, which have not been studied in detailed yet. More 
experimental input is needed to study these multi-body decay modes. 

Study of more than three-body B decay modes with Xci and Xc 2 in final state 
not only help in understanding B meson decays but also provide portal to search for 
charmonium/charmonium-like exotic states in one of the intermediate hnal states. 
For example, looking at the Xcit^~^'^~ invariant mass spectrum in i? —)■ Xci'^~^'^~K 
decays, one can search for Xci(2P) and/or X(3872). 

Using the Xcj J/'4^1 modes, we report on the inclusive branching fraction (B) 
of B —)■ XcjX decays. To further understand xa and Xc2 production in B decays, 
we reconstruct the following exclusive B decays: B^ —)■ Xcj'^~K^, 

B^ Xcjt^^K+, B+ and 5° [7]. 


2 Data sample and event selection 

We use a data sample of 772 x 10® BB events collected with the Belle detector [ 8 ] 
at the KEKB asymmetric-energy e’''e“ collider [9] operating at the T(4S') resonance. 
All results presented here are preliminary. 

The meson is reconstructed via its decays to (£ = e or p). To reduce the 
radiative tail in the e“*'e“ mode, the four-momenta of all photons within 50 mrad with 
respect to the original direction of the or e~ tracks are included in the invariant 
mass calculation, hereinafter denoted as Me+e-( 7 ). The reconstructed invariant mass 
of the J/4j candidates is required to satisfy 2.95 GeV/c^ < Me+e-(-^) < 3.13 GeV/c^ 
or 3.03 GeV/c^ < M^+^- < 3.13 GeV/c^. For the selected J/i^ candidates, a vertex- 
constrained fit is applied and then a mass-constrained fit is performed in order to 
improve the momentum resolution. 

The Xci and Xc2 candidates are reconstructed by combining J/'ip candidates with 
a photon having energy {E^) larger than 100 MeV in the laboratory frame. To reduce 
the combinatorial background coming from vr® —)■ 77 , we use a likelihood function that 
distinguishes an isolated photon from vr® decays using the photon pair invariant mass, 
photon laboratory energy and polar angle HD]. We reject both 7 ’s in the pair if the 
TT® likelihood probability is larger than 0.3 (0.8) for an inclusive study {B —)■ XcjKtt, 
B XcjEtiti decay modes). 
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3 Inclusive B decays to XcJ 

To identify the signal, we use the J/V'y invariant mass Mj/^^ distribution and extract 
the signal yield from a binned maximum likelihood ht. A double-sided Crystal Ball 
function is used to model the signal shapes of i? —?• XciX and B Xc 2 ^- Figured] 
(a) shows the ht to the distribution for B —)■ XciX and B —)■ Xc 2 ^ decays. 

The reconstruction efhciencies for the inclusive B —)■ Xci^ and B —)■ Xc 2 ^ decays 
are estimated to be 24.2% and 25.9%, respectively. Uncertainty on the efficiency is 
estimated to be 4.0%. 
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Figure 1: (a) distribution of the B — Xcj{^ decays in data. 

The curves show the signal (cyan dash-dotted for Xci and red dashed for Xc 2 ) and the 
background component (green dash-double-dotted for combinatorial) as well as the 
overall ht (blue solid). The lower plot shows the pull of the residuals with respect 
to the ht. (b) Plots showing DB{B —Xci^) and DB{B Xc 2 ^) in each bin of 
for B —> Xcj{^ These plots are without feed-down subtraction. The 

uncertainty shown in these plots are statistical only. 


After subtracting the feed-down contribution, we get the direct branching 
fractions B{B Xci^) and B{B Xc 2 X) to be (3.03 ± 0.05 ± 0.25) x 10“^ and 
(0.70 ± 0.06 ± 0.10) X 10“^, respectively. First (second) error is statistical (system¬ 
atic). Here, the systematic uncertainty dominates the measured branching fractions. 

Figured] (b) shows the obtained distribution of the DB in bins of Suppression 
of two-body decay of Xc 2 is visible in the p^j distribution. Most of the Xc 2 production 
comes from decays of three-body or more. 
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4 Exclusive reconstruction 


The reconstructed invariant mass of the Xci (Xc 2 ) is required to satisfy 3.467 GeV/c^ < 
< 3.535 GeV/c^ (3.535 GeV/c^ < < 3.611 GeV/c^). A mass-constrained 

ht is applied to the selected Xci and Xc 2 candidates. 

XcJ candidates are combined with charged kaon and pion candidate tracks to re¬ 
construct B meson. To identify the B meson, two kinematical variables are used: 
beam-constrained mass (Mbc) and energy difference (Ai7). The Mbc is dehned as 

\J^hea.ra ~ and the AE is defined as Ei — A'beam, where i^beam is the beam 

energy in the GM frame and Pi {Ei) is the momentum (energy) of the i-th daughter 
particle in the GM frame and the summation is over all final states used for reconstruc¬ 
tion. We reject candidates having Mbc less than 5.27 GeV/c^ or |Ai7| > 120 MeV. 
We extract the signal yield from an unbinned extended maximum likelihood (UML) 
ht to the AE variable. Fig. |2] and |3] shows the ht to the AE distribution for the decay 
modes of interest. 




Figure 2: AE distribution for (a) B^ Xcit^ (b) B^ —)■ Xc 2 '^ , (c) B^ 

(d) B+ Xc27r+iF5, (e) B+ XciTr°K+, and (f) B+ -)> Xc 2 vr°iF+ decay 
modes. The curves show the signal (red dashed), peaking background (magenta dash- 
dotted) and the background component (green dotted for combinatorial) as well as 
the overall ht (blue solid). 


To understand the production mechanism of intermediate states, we look at the 
background subtracted Mk-k-k, and M^^^ distributions for the 

decay mode of interest. We perform a UML ht to the AE distribution and use the 
sVlot formalism m to project signal events in the distribution of interest. 
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Figure 3: AE distributions for (a) —)■ (b) B'^ —)■ Xc 2 T^^t^ K^, 

(c) 5° -)> Xci7r+7r"iF^, (d) 5° -)> Xc27r+7r"iF^, (e) 5° -)> Xcivr°7r"7F+ and (f) 
—)■ Xc27r°7r“iF+ decay modes. The curves show the signal (red dashed), peak¬ 
ing background (magenta dash-dotted) and and the background component (green 
dotted for combinatorial) as well as the overall £t (blue solid). 


The K*{892) is found to be a major contribution in the B —)■ Xci^K decay modes 
as seen from Fig. 0] (a), (e) and (i), while in i? —)■ Xc 2 '^K decays the iF*(892) com¬ 
ponent is less prominent. Our study suggests that the B —)■ Xc2K*{892) mechanism 
does not dominate the B Xc 2 '^^ decay, which is in marked contrast to the Xci 
case. From this study one can say that the production mechanism of the Xc 2 from B 
mesons is different in three-body decays for the B —)■ Xcjt^K case. 

Background subtracted s'Plot distribution of Mk-k-w, Mk+tt-, and 

are shown in Figs. |5]and[6]for the B^ —)> Xcj'^^'^~K^ decay mode. No narrow 
resonance is seen in the and distributions with the current statistics. 

In our search for the A(3872) and/or Xcii!^P) decaying into we didn’t 

hnd any signal and provided upper limit (@ 90% C.L.) as: 

• B(B± ^ A(3872)iF±) x B(A(3872) ^ Xcivr+vr-) < 1.4 x lO”® 

. B{B+ ^ x,i(2P)iF+) X B(xei(2P) ^ Xei(lP)vr+7r-) < 1.1 x lO’^ 

5 Summary 

Belle presented the preliminary results at this conference. We measured the direct 
branching fractions B{B —)■ XaX) and B{B —)■ Xc 2 X) to be (3.03 ±0.05 ±0.25) x 10“^ 
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Figure 4: Background subtracted s'Plot Mk-k and distributions for (a and b) 

—)■ XciTT^iF’*', (c and d) —)■ Xc 2 't^~K^^ (e and f) B~ —)■ Xci^r^iF^, (g and h) 

B~ —)■ Xc 2 '^~Kg and (i and j) B~ Xci^^^K^ decay modes. 
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Figure 5: Background subtracted s'P^ot (a) (b) (c) and 

(d) M^^ 27 r± distributions for 5+ —)> Xcjt^^t ^~decay modes. 
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Table 1: Summary of the results. Signal yield (T) from the £t, signihcance (S) 
with systematics included, corrected efficiency (e) and measured jB. For jB, the hrst 
(second) error is statistical (systematic). Here, in the neutral B decay case, the 
—)■ TT+TT” branching fraction is included in the efficiency (e) but the factor of 2 (for 
—)■ or iF°) is taken into account separately. TZq is the ratio of B{B —)■ Xc 2 ^) 

to B{B —)■ Xci^), where X is the same set of particles accompanying the Xci (Xc 2 ) in 
the hnal states. 


Decay 

Yield (Y) 

5(a) 

6(%) 

B (10-4) 

71)3 

B^ ^ Xcj7r-K+ 




0.14 ±0.02 

Xcl 

2774 ± 66 

66.7 

17.9 

4.97 ±0.12 ±0.28 


Xc2 

206 ± 25 

8.7 

16.2 

0.72 ±0.09 ±0.05 






0.20 ±0.04 

Xcl 

770 ± 35 

33.7 

8.6 

5.75 ±0.26 ±0.32 


Xc2 

76.4 ± 14.7 

4.6 

7.5 

1.16 ±0.22 ±0.12 


B+ ^ Xcj7r^K+ 




< 0.21 

Xcl 

803 ± 70 

15.6 

7.8 

3.29 ±0.29 ±0.19 


Xc2 

17.5 ±28.4 

0.4 

7.0 

< 0.62 


B"^ XcjTT+vr 




0.36 ±0.05 

Xcl 

1502 ± 70 

19.2 

12.8 

3.74 ±0.18 ±0.24 


Xc2 

269 ± 34 

8.4 

11.4 

1.34 ±0.17 ±0.09 


B^ -)■ XcjTT+vr 




< 0.61 

Xcl 

268 ± 30 

7.1 

5.4 

3.16 ±0.35 ±0.32 


Xc2 

37.8 ± 14.2 

1.8 

4.8 

< 1.70 


B° Xcjvr°7r K+ 




< 0.25 

Xcl 

545 ± 81 

6.5 

5.0 

3.52 ±0.52 ±0.24 


Xc2 

-76.7 ±42.0 

- 

4.3 

< 0.74 
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Figure 6: Background subtracted s'Plot (a and d) Mk+tt+-k-, (b and e) Mx+n- and 
(c and f) M^+^+ distributions for i?+ decay (upper) and —)■ 

Xc 2 '^~^'^~K~^ decay (lower), respectively. 

and (0.70 ± 0.06 ± 0.10) x 10“^, respectively. We observe the Xc 2 '^~K~^ decay 

mode for the hrst time, with 206 ± 25 signal events having a 8.7a signihcance, along 
with evidence for the B~^ Xc 2 '^^K^ decay mode, with 76±15 signal events and 

a signihcance of 4.6cr. In four-body decays, we observe the B^ —)■ Xci7r'''7r“iF’'', 
B^ -X- XciT^^T^~, B^ —)■ Xci7r''“7r“iF5, and B^ —)■ Xci7r°7r“iF’'' decay modes for the 
hrst time and report on measurements of their branching fractions. We hnd that Xc 2 
production in comparison with Xci increases with a higher number of multi-body B 
decays. We observe that the Xc 2 is more often accompanied by higher K* resonances 
as opposed the Xci which is dominantly produced with lower K* resonance. Inclusive 
and exclusive study of B decays having Xc 2 in the hnal state suggests suppression of 
two-body B decay due to suppression of a tensor, while multi-body B decays into 
Xc 2 are allowed. In our search for X(3872) —)■ and Xci(27^), we determine 

an U.L. on the product of branching fractions. The negative result for our searches 
is compatible with the interpretation of X(3872) as an admixture state of a D^D*^ 
molecule and a Xci{‘^P) charmonium state. 
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